Role of hyperfine interaction for cavity-mediated coupling between spin qubits 
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We consider two qubits interacting by means of an optical cavity, where each qubit is represented 
by a single electron spin confined to a quantum dot. It is known that electron spins in III-V 
semiconductor quantum dots are affected by the decoherence due to the hyperfine interaction with 
nuclear spins. Here we show that the interaction between two qubits is influenced by the Overhauser 
field as well. Starting from an unpolarizied nuclear ensemble, we investigate the dependance of 
the fidelities for two-qubit gates on the Overhauser field. We include the hyperfine interaction 
perturbatively to second order in our analytical results, and to arbitrary precision numerically. 



I. INTRODUCTION 

Substantial progress has been made towards the imple- 
mentation of coherent interfaces between an electron spin 
in a quantum dot and photons [1J44] . The coupling of sin- 
gle spins to photons is a promising mechanism for imple- 
menting quantum information processing schemes based 
on quantum dots, where a qubit (quantum bit of infor- 
mation) is represented by electron spin degrees of free- 
dom [5], The main source of electron spin decoherence 
in semiconductor quantum dots is considered to be the 
hyperfine interaction to the nuclear spins of the lattice 
atoms O [7J . The dephasing time due to the interaction 
with an unprepared ensemble of nuclear spins is about 
10 ns [7J which is much shorter than the single qubit 
operation times by means of electronic control (« 100 ns 
[9]). There are some ways to deal with this problem: to 
prolong the decoherence time using advanced spin-echo 
techniques |10j or dynamical nuclear polarization [11] , to 
use nuclear spin "poor" materials such as graphene |12j 
or nitrogen- vacancy centers in diamond |13j or to shorten 
the qubit manipulation time [14 . The optical driving of 
electron spins has this last advantage and offers mecha- 
nisms for single qubit operation time around 30 ps [2]. 
Besides single qubit manipulations, two qubit operations 
are required for implementing quantum algorithms. One 
of the schemes for optical exchange interaction between 
two spins is based on virtual Raman transitions between 
valence band heavy holes and conduction band electrons 
[15j . Thereby, a set of quantum dots is coupled to a 
high finesse cavity. Each quantum dot is doped with a 
single electron and can be controlled by linearly polar- 
ized laser light. The laser polarization is perpendicular 
to the growth direction of the quantum dots and perpen- 
dicular to the cavity field polarization as well, so that 
laser light is acting on the quantum dots without driv- 
ing the cavity. A magnetic field is applied in the Voigt 
configuration, i.e., perpendicular to the growth direction 
and parallel to the cavity field polarization, so that single 
photons emited by quantum dots escape into the cavity 
and can be detected afterwards. The cavity frequency is 
chosen such that it is slightly detuned from the transi- 
tion from the electron spin down state to the trion state, 
while the laser frequency is detuned from the transition 
from the spin up state to the trion state. If these detun- 




Figure 1. Cavity mode mediated spin-spin interaction. If the 
difference of the detunings A for two quantum dots coincides, 
then they can exchange virtual cavity photons of energy w c . 
Here, uj c and io\ are the frequencies of the cavity mode and the 
laser acting on quantum dot i, and A denotes the difference 
between the detunings, i.e. the detuning from the two-photon 
resonance on a given quantum dot. 



ings are the same then the A system is in a two-photon 
resonance that allows for manipulating a single electron 
spin (Figs. [I] and [2]). If the detunings in a quantum dot 
are different, but the difference of the detunings (the de- 
tuning from the two-photon resonance) match for a pair 
of quantum dots, then the electron spins in these quan- 
tum dots can resonantly exchange energy through virtual 
cavity photons (Fig. [IJ. In this way, virtual photons 
can be used for performing coherent interaction between 
arbitrary qubits. Two-qubit operations can potentially 
be carried out on a subnanosecond timescale |T5j, there- 
fore the decoherence can be expected to originate mainly 
from nuclear spins. In this paper, we calculate the fideli- 
ties of such two-qubit operations in the presence of the 
hyperfine interaction. We will study in particular par- 
tial SWAP operations and the controlled-not (CNOT) 
gate. For this purpose, we first calculate the fidelity of 
the quantum gates of interest analytically to second or- 
der of the Overhauser field and subsequently we perform 
an exact numerical calculation for comparison. 



II. MODEL 

We start from the microscopic model describing a 
quantum dot embedded in a cavity and irradiated by a 
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laser field, 



with the anti-Hcrmitian operator 
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where ui a and e CT are the energy and annihilation opera- 
tors for the electronic states in the conduction band with 
spin a = ±1/2 =f, -J, and in the valence band with total 
angular momentum a = ±3/2 (we assume that the heavy 
hole subband is sufficiently split from the light hole band 
to allow for pure heavy hole excitation). The quantiza- 
tion axis for the electrons is set by the external magnetic 
field along the x direction, while the quantization axis 
in the valence band is given by the structure and as- 
sumed to be in z direction. Furthermore, for a linearly 
polarized cavity mode (w c , a c ) along the x direction and 
linearly polarized laser mode (wj,,az,) in the y direction, 
we obtain a radiative coupling to the linear combination 
e v = (e_ 3 /2 — e3/2)/v / 2 in the valence band. The cou- 
pling strengths to the cavity and laser modes are denoted 
with g c and <7l, respectively. The hyperfine interaction 
of a conduction band electron in the quantum dot with 
the surrounding nuclear spins can be written as 
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where I k is the operator belonging to the fc-th nuclear 
spin in contact with the electron and Ak denotes the 
corresponding hyperfine coupling constant. The aver- 
age hyperfine coupling constant for GaAs is 90 fj,eV 
|16j . The electron spin operator is given by S = 



2 L-ia.a 1 - 



where cr = (<r x , a y ,a z ) is the vec- 



tor consisting of Pauli matrices. Here, we neglect the 
dipolar hyperfine coupling of the valence band states 
since it is typically smaller. In the presence of a magnetic 
field along the x direction exceeding the nuclear field of 
typically ~ 10 mT, we can neglect the flip-flop terms and 
write Hhf rj S x Ylk=i Aklx = S x h, where h denotes the 
Overhauser (nuclear) field operator in x direction. We do 
not take the nuclear Zeeman terms into account, because 
they are considerably smaller than the thermal energy 
and the electron Zeeman energy. We can then combine 
the hyperfine Hamiltonian with the first term in Eq. ([!]) 
by using oj a = (g e fiBB + h)a/2 for a =t,i= ±1, where B 
is the magnetic field applied along x and g e is the effec- 
tive electron g-factor. We can now replace the operator 
h with one of its eigenvalues and perform an average over 
h later. This allows us to follow the steps performed in 
Ref. |15) before taking the average over nuclear configu- 
rations. In GaAs quantum dots the number N of nuclear 
spins is large, typically between 10 5 and 10 6 , and there- 
fore the Overhauser field follows a Gaussian distribution 
around mean value and with variance a ~ A/yN [15] . 

In order to eliminate the valence band states, we per- 
form a Schrieffer- Wolff transformation H e g = e~ s He s 
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where we use W3/2 ~ — ^v- By setting the va- 

lence band occupation number to e\e v ps 1 and using 
that the laser field is in a coherent state with g^o-L ~ 
VLl exp(— iu)i,t), we obtain to second order in the cavity 
and laser couplings g c and g^, 
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where the index i refers to the respective quantities in 
dot i — 1,2, and dij = i,j =t>4- The first term 

inside the sum in Eq. Q is the Zeeman splitting due to 
the external and nuclear fields, the second term describes 
the effective coupling of the cavity mode to the quantum 
dot electron spins with strength 
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are the detunings of the cavity and laser fields (s. Fig. 
PI), and the last two terms in Eq. (|| can be interpreted 
as the Lamb and Stark shifts of the cavity and quantum 
dot levels, respectively. 

A second Schrieffer- Wolff transform can be used to also 
eliminate the cavity mode, which leads to the effective 
photon-mediated interaction between two spins i and j 
in the interaction picture with 
H o = J2t(9eHBB + hi)a z n , 
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Ai{hi) = Ai{hi) - Ai{hi) 
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represents the coupling strength between two electron 
spins. The time-independent interaction Hamiltonian 
Eq. ^ strictly applies to the two-photon resonance 
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Figure 2. Energy level scheme for a quantum dot filled with 
a single electron and coupled to a cavity mode. The Zeeman- 
split single-electron states can be excited to a trion (negatively 
charged exciton) state \X~) by coupling to the cavity or laser 
field. Both cavity and laser field frequencies are detuned by 
A c and Al from resonance, and the combined system is de- 
tuned from its two- photon resonance by A = A c — Al- The 
Overhauser shift caused by the hyperfine coupling to the nu- 
clear spins leads to a fluctuating detuning, thus reducing the 
fidelity of the optically generated quantum gates. 



Ai = Aj = A in the absence of nuclear spins. By going 
over into the rotating frame with nuclear spins by R — 
exp (it(hi<7^ + hjoLS)) we obtain a time-independent ef- 
fective interaction Hamiltonian with Heisenberg trans- 
verse coupling type between two electron spins: 



The unitary time evolution operator of the interaction 
between two electron spins is 



U {<f)) = exp 



HlLdt 
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Figure 3. Effect of the Overhauser fields hi 
time evolution of the two-electron spin state: change of the 
interaction phase and precessing of the rotation axis in the 
x-y plane with (hi — hj)t 



two-electron spin state affected by nuclear spins coincides 
with an unaffected one. 

The two-qubit CNOT operation can be implemented 
as a sequence of single spin rotations combined with the 
unitary time evolution operator [15 : 
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with cf> — J cjij (t) dt. The time evolution operator in orig- 
inal frame is then 
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In the subspace {| f4-)) I it)} the operator U(4>) acts 
as a rotation exp(i^xr • h) with n = (— cos{(ft-i — 
hj)t},sin{(hi — hj)t},0). Thus the changed time evo- 
lution of the two-electron spin state due to hyperfine 
interaction can be interpreted as a modified rotation 
in the mentioned subspace (s. Fig. [3]). There are 
two distinct effects due to the nuclear spins. First, the 
interaction phase becomes Overhauser field dependent 
4> = 4>{hi,hj) and, secondly, the rotation axis starts to 
precess in the x-y plane with (hi — hj)t. The second 
effect has the maximal contribution when 0(0) = wf , 
where 0(0) = gij(hi — 0,hj = 0)t = gij(0)t and where 
m is an integer. When 0(0) = mir, the trajectory of the 



where hi = (1, 1, — l)/v3, hj = ( — 1, 1, l)/y/3. 

We calculate the fidelity of the generated unitary gate 
in the presence of the Overhauser field with respect to the 
ideal unitary gate (without Overhauser field) in order to 
quantify the effect of the nuclear spins. This fidelity re- 
flects the difference in final state after the gate operation, 
averaged over pure input quantum states |17j . The aver- 
age fidelity for operators acting in 4-dimensional Hilbert 
space is given by [18] 



F(d(o),d(h)) = 



|Tr(Ot(0)Q(fe))|; 
20 



(15) 



where 0(h) and O(0) denote the unitary operator with 
and without Overhauser field dependence. Next, the av- 
erage fidelity Eq. ( p~5| ) must be averaged over the nuclear 
spin field distributions in the quantum dots i and j to 
extract the Overhauser field dependence of the fidelities 
for the time evolution and the CNOT operator. 
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III. FIDELITIES IN SECOND ORDER 
HYPERFINE INTERACTION 

Starting from the microscopic model Eq. and per- 
forming two subsequent Schrieffer- Wolff transformations 
[IS] while including the hyperfine interaction, we find the 
effective spin-spin interaction Hamiltonian as described 
above. Using Eqs. Q, and ([To|), we find 

that the effective spin-spin coupling depends on the Over- 
hauser fields hi and hj in the quantum dots, 
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We assume that the fluctuations of the Overhauser field 
around its zero mean value are small, since h ~ A ~ 10~ 5 
eV. Therefore we can investigate analytically the fideli- 
ties to the second order of the Overhauser field. Since 
the interaction Hamiltonian ([8| is nuclear field depen- 
dent partly through the coupling coefficient , we need 
to expand it with respect to the Overhauser field in both 
coupled quantum dots, 
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The trace of the product of the perfect time evolution 
operator and the one with included hyperfine interac- 
tion is Tr[E/t(o) U(h)] = 2{1 +cos[^(0)] cos(0) + cos[(^ - 
hj)t] sin[</>(0)] sin(0)}. We obtain for the average fidelity 
between those two operators to the second order of the 
Overhauser fields hi and hj, 

F(U)=l- ^(OXD^+D^ 2 ) 

-\ ^\m\'W-Ah l hj? + 0(h 3 ), (18) 
9ij\y) 

where D ; = (A/A l c A l L - l/A)/2, I = are the first- 



order coefficients in the Overhauser field in Eq. (17 1, and 



the time was expressed through the interaction phase 
0(O) = ft 3 -(O)t. 



Averaging of the fidelity ( 18 1 over both Overhauser 
fields, we find 



{F(U)) hi , hj 
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sin 2 [0(O)]cr 2 + 0(h 3 ). (19) 



Assuming the parameters Al = 4.5 meV, A c = 5 meV, 
g c = 0.5 meV, fl l L = £l J L = 1 meV and A = 0.5 meV, such 
that D = D, = D., = -0.99 x 10 3 eV _1 and g M = 0.02 



meV, the average fidelity given by Eq. ( 19 ) behaves as 
function of <p(0) as shown in Fig. [4] The value of (F(U))h 
by <j){0) = tt/4 is equal to 0.99995. 

As we can see from the formula ( [19] ), the two effects 
induced by nuclear spins decrease the fidelity with cor- 
responding strengths: the changing of the electron spin- 
spin coupling strength and resonance conditions oc D 2 
and the precession of the rotation axis of the two-electron 
spin state oc l/g 2 , (0). These are determined in Eqs. (|5|), 

uch 



(16), and (17) and depend on tunable parameters sucr 
as A and A^ (we can express A c = A^ + A). The cav- 
ity detuning can be tuned accordingly, e. g. in photonic 
two-dimentional slab microcavities using Xe condensa- 
tion [TSJ. The contributions to fidelity reduction from 
terms cx D 2 and cx l/gfj(0) are not of the same order (s. 
Fig.|5]b), D 2 < 0.1/5 2 (0). But we still keep the terms 
oc DMn our analytical calculations, because they give the 



upper boundary of the fidelity, when sin[</>(0)] = 0. By 
changing the laser detuning A^ and the detuning A we 
can adjust the values D and cjij and thus tune the fidelity 
itself. 

The dependence of the fidelity (for <fi(0) = ~) in respect 
to A and A^ is shown in Fig. [5]ja). It indicates that fi- 
delity can be improved considerably by changing the de- 
tunings. However by choosing better parameters, some 
assumptions of the applied theory should be conserved, 
so that gij <C A and Al(A c ) » g c . One can tune the 
intensity of applied control laser and rotate the electron 
spins faster, so that they interact faster as well and "feel" 
the dephasing in a lesser extent. However the coupling 
strength (jij(t) depends quadratically on laser intensity 
and just slight enhancement of Rabi frequency Ql will 
make it of the same order as A. Therefore, to improve the 
fidelity and to keep these relations true one can use e. g. 
A = 0.3 meV and A^ = 4 meV. With this choise, the av- 
erage fidelity for the gate U(j ) is increased to 0.999991. 
The fidelity decays considerably slower in this case as 
well (Fig. [4]). The formula Eq. (15) can be applied for 



calculating the average fidelity of the CNOT operation. 
If we neglect the imperfections of the single qubit opera- 
tions in the CNOT operator sequence and focus only on 
the decoherence due to the interaction part, then we find 
for the trace distance between a perfect CNOT operator 
and CNOT operator with nuclear spins, 

Tr [^CNOT #CNOt(M] = 
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The CNOT average fidelity in second order in the Over- 
hauser field is then 
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Figure 4. Fidelity F of the unitary time evolution operator 
U generated by the XY interaction in the presence of nuclear 
spins as a function of the interaction phase <j>(0) oc t, where t 
denotes the interaction time. Both lines show the analytical 
result calculated in second order perturbation theory in the 
nuclear field. The fidelity can be increased by tuning the 
parameters A and Al- For smaller A and At, the fidelity 
increases due to faster gate operation. The other parameters 
are held fixed and are given in the text. 
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Figure 5. a) The average fidelity F of the unitary time evolu- 
tion operator U for <^>(0) = n/4 calculated in second order of 
the nuclear field as function of the laser detuning Al and two- 
photon detuning A. An adjustment of these parameters can 
lead to an increase of the fidelity or to its reduction. However, 
the choice of the optimal parameters is limited by the condi- 
tions of the applied formalism. Further details are given in 
the text b) The ratio between two contributing mechanism to 
reduction of the fidelity D 2 (A, Al) ■ glj{A, Al)- It increases 
for smaller detunings and consequently for improved fidelities. 
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Figure 6. The analytical and numerical results for fidelity of n 
subsequent CNOT operations in the presence of nuclear spins 
for different detunings. The analytically calculated values are 
represented by circles and the numerically calculated ones by 
stars. 



We find that F CNOT = 0.9999 by using A = 0.5 meV 
and A L = 4.5 meV, and F CNOT = 0.999983 with A = 
0.3 meV and Al = 4 meV, which is smaller than the 
fidelity of a single U(j ) operation. 

For the average fidelity for further applied CNOT op- 
erations, we find 
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CNOT 
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The fidelity F CNOT (n) averaged over the nuclear field 
(Fig. [6]) is reduced after each application and can be 
increased by a better choice of the detuning parame- 
ters. Thus for the last set of parameters only after seven 
CNOTs, the fidelity drops below typical error correction 
thresholds (< 0.9999) gU]. 



IV. NUMERICALLY EXACT FIDELITIES 

For large values of the parameter <fi(0) we evaluate the 
averaging of the fidelities over nuclear fields without ex- 
panding |Tr([/(0) t ?7(/i))| 2 . The results of the numerical 
averaging are shown in the Fig. [7j The comparison to 
the analytical result in Fig. ^&) reveals that our analyt- 
ical model can predict the behavior of the fidelity only 
up to a certain value of </>(0). From Fig. [7ja) one can 
see that around </>(0) ~ 12 the numerical and analytical 
data begin to differ. For large 4>(0), the analytical data 
even diverges into negative values. The numerical data 
on the other hand converges to 1/2. The numerical eval- 
uation of the fidelity for different sets of detunings [Fig. [7] 
(b)] confirms that the tuning of the transition parame- 
ters can improve the fidelity significantly. The numerical 
curve oscillates with the same period 27T as the analyti- 
cal one. The maxima correspond to the situation, when 
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Figure 7. Numerically calculated average fidelity of the time 
evolution operator for large 0(0). (a) The comparison be- 
tween analytical and numerical calculations for = 4.5 meV 
and A = 0.3 meV. The analytically obtained fidelity diverges 
for large 0(0) while the numerical fidelity saturates to 0.5. (b) 
Comparison between two numerically calculated fidelities for 
the gate U((f>) with different detunings. 



the fidelity is decreased only due to the modified inter- 
action phase between the electron spin and there is no 
effect of the changed rotation axes of the two-electron 
spin state. The maxima correspond to the decay of the 
fidelity, where both Overhauser field induced effects are 
contributing and the effect of the precession of the rota- 
tion axis is maximal. 

The fidelity curve suggests that the behavior of its 
lower boundary could be described by the Gaussian 

e -0(o) _ The fidelity envelope to the second order 
of 0(0) is obtained from Eq. (19) by setting sin(0(O)) 2 
= 1, We fit this envelope to the decay function |[1 + 
exp(— 0(O) 2 /t 2 )], where r is given through parameters 
a, D and c/ij from Eq. ( 19 ) and is defined as 



9iM 
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r(A,A L ). 
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The Gaussian function gives a good fit to the minimum 
points of the fidelity for small (f>(0), but decays faster 
for larger (f)(0). The tail of the fidelity can be fitted by 



adding a third-order term to the argument of the expo- 
nential function, j[l+exp(— 0(O) 2 /r 2 +c0(O) 3 )] withe as 
a fitting parameter. The fits are shown in Fig. ^ and fit 
with T\ and c\ corresponds to the curve with parameters 
Al = 4.5 meV and with t 2 and c 2 to numerical curve with 
another set of parameters. The fits yield c\ = 9.6 x 10~ 7 
and c 2 = 5 x 10~ 8 . 

The time dependance of the fidelity can be found by 
changing from interaction phase 0(0) to time by assuming 
0(0) = (jij(0)t. Then for short time dependance of the 
lower boundary of the fidelity is described by the function 
Q.5[l + exp(-t 2 /T 2 )], where 
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Since from Fig. (p| we know that D 2 g 2 ) (0) < 1, we can 
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neglect it here. It follows that T = 

where N is the number of the nuclear spins and A is 
hyperfine interaction constant. This implies that the de- 
cay of the fidelity does not dependend on the interaction 
strength of the electron spins and is only affected by the 
nuclear field distributions. By tuning the transition pa- 
rameters we increase the interaction strength between the 
coupled electron spins to some reasonable value. Thus 
the spins acquire a larger interaction phase within the 
same time. But the decay time of the time evolution op- 
erator fidelity is the same and for GaAs (N = 10 5 , A = 
90 fieV) T = 2 ns. The T 2 decay time of a electron spin 
coupled to the nuclear bath is given by T 2 = 2hy ^ [2Tj 
and in GaAs is ~ 5 ns. Since both interacting spins are 
affected by the decoherence, the decay of the common in- 
teraction phase is faster than of a single one (T w 0.4 T 2 ). 

The numerically exact CNOT fidelities are presented in 
Fig. [6j The second order approximation in Overhauser 
field gives a good agreement to the numerical results. 
Since the CNOT operator is defined for a fixed interaction 
phase, there is no oscillation behavior and the fidelity 
decreases for each further applied CNOT gates. 



V. CONCLUSIONS 

We calculated average fidelities for unitary time evolu- 
tion operator and CNOT operation of two quantum dots 
interacting through an off-resonant cavity mode. We ob- 
tained results to second order of hyperfine interaction 
analytically and to an arbitrary order by averaging the 
fidelities over the Overhauser field distributions numeri- 
cally. Both approaches are in good agreement for small 
interaction time durations. If only the longitudinal com- 
ponent of the Overhauser field in two interacting quan- 
tum dots is considered for the hyperfine interaction, then 
the time evolution of the two-electron spin state changes 
in two ways (Fig. [3]), and both effects contribute differ- 
ently to the reduction of the fidelity. A prolongation of 
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the decoherence time of the two-electron state is possi- 
ble by narrowing the distribution of the Overhauser field 
22J [Eq. (19), (21)]. The precession of the rotion axis of 
the two-electron state on the Bloch sphere can be min- 
imized by preparing the the nuclear spin states in each 
interacting quantum dot in the same certain state and 
minimizing the fluctuations of the Overhauser fields [33J . 
Since we neglected the hyperfine spin-flip terms in our 



calculations, the obtained results are only valid for the 
Zeeman fields oc 1 T. 
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